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Quantum Logic with Interacting Bosons in ID 
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We present a scheme for implementing high-fidelity quantum logic gates using the quantum walk of 
a few interacting bosons on a one-dimensional lattice. The gate operation is carried out by a single 
compact lattice described by a one-dimensional Bose-Hubbard model with only nearest-neighbor 
hopping and on-site interactions. We find high-fidelity deterministic logic operations for a gate set 
(including the CNOT gate) that is universal for quantum information processing. We discuss the 
applicability of this scheme in light of recent developments in controlling and monitoring cold-atoms 
in optical lattices, as well as an implementation with realistic nonlinear quantum photonic devices. 


Introduction — Quantum walks (QWs) are unitary pro¬ 
cesses that describe the quantum-mechanical analogue 
of the classical random walk process m- Since their 
conception, there has been a broad interest in their pos¬ 
sible use for quantum information processing fa in. Two 
mathematical models for QWs have been developed: the 
discrete-time QW |T| , in which the particle takes discrete 
steps in a direction given by a dynamic internal degree 
of freedom (a coin), and the continuous-time QW [5] in 
which the dynamics are described by Hamiltonian evo¬ 
lution on a lattice in the tight-binding representation. 
Here, we consider the continuous-time quantum walk on 
a one-dimensional lattice. 

Experimentally, QWs have been implemented with pho¬ 
tons EHm, trapped ions maun, and ultra-cold atoms 
GMZI, among other platforms. Specifically in the field 
of ultra-cold atoms dim, the degree of experimental 
control is remarkable: it is possible to prepare an ini¬ 
tial state with single-site and single-particle resolution, 
to create arbitrary one- or two-dimensional lattice po¬ 
tentials, to determine the interaction between the parti¬ 
cles, and to directly monitor in real space the evolving 
many-body distribution. 

The ability to control and monitor quantum particles 
with such precision offers an interesting route to the 
implementation of quantum information processing and 
quantum computation schemes. Universal quantum com¬ 
putation has been theoretically shown possible using 
QWs with interacting particles on certain non-trivial two- 
dimensional lattices [HI ESI and on one-dimensional lat¬ 
tices with a large number of degrees of freedom at each 
lattice site |21If25] . However, an implementation using 
quantum particles hopping on a simple one-dimensional 
lattice, without any additional degrees of freedom, is not 
known. Such a geometry would greatly simplify the ex¬ 
perimental implementation, bringing it into the realm of 
recently reported experimental techniques d- Further¬ 
more, one-dimensional implementations offer other im¬ 
portant practical advantages, for example freeing the sec¬ 
ond spatial dimension for important tasks such as error 
correction or connecting remote qubits. Other important 


Figure 1: (color online). Illustration of one-dimensional quan¬ 
tum walk based quantum gates, (a) A qubit in the state |0) 
with dual-rail encoding, (b) Implementation of a single-qubit 
gate, (c) Schematic of a two-qubit system on a lattice. 

tasks such as process tomography could still be performed 
in ID (see Supplementary Information). 

In this work, we show how it is possible to use multi¬ 
particle continuous-time quantum walks in a simple ge¬ 
ometry — a one-dimensional lattice with only nearest- 
neighbor hopping and on-site interactions — as a com¬ 
pact platform for implementing quantum logic. We 
demonstrate our approach by detailing a set of lattice 
potentials that yield, with high fidelity, a universal set 
of quantum gates with only two sites per qubit. More¬ 
over, the required lattice potential for each gate is time- 
invariant, simplifying the experimental implementation 
and possibly reducing the total operation time. Thus, 
high-fidelity gates can be constructed with lower prob¬ 
abilities of qubit-loss errors (i.e. lost particles during 
the computation) that necessitate computationally costly 
error-correction procedures. While we focus on interact¬ 
ing ultra-cold bosonic atoms trapped in an optical lattice 
— a physical system in which our results can be imple¬ 
mented with existing experimental techniques m — we 
also discuss the possibility of extending our analysis to 
nonlinear quantum photonic systems. 

The dynamics of bosonic particles on a lattice is 
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described by the time-independent many-body Bose- 
Hubbard Hamiltonian 


H=Y J En 
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where E m is the on-site energy of site m, aj Tl \a m is 
the creation\annihilation operator for a boson in site m, 
n m = alndm is the number operator, < 0 is the 

tunneling rate between nearest neighbors, and T is the 
on-site interaction energy that arises when two or more 
bosons occupy the same site. The unitary transformation 
describing the evolution of multiple quantum particles 
propagating on the lattice is given by e ~ lHt , where t is 
the propagation time. The quantum logic gates discussed 
here will be implemented by evolving under this Hamil¬ 
tonian (with a suitable choice of parameters) for some 
predefined time tfi na i which we take to be tfi na i = 1. 


Defining qubits on a lattice — The basic element of in¬ 
terest for quantum gates of the type discussed here is the 
quantum bit, or qubit. The continuous-time quantum 
walk, however, is described by the evolution of quantum 
particles on a lattice according to the Hamiltonian de¬ 
scribed in Eq. |T] To define our qubits on the lattice, 
we use a spatial encoding where a qubit is physically 
implemented by a single boson in a pair of neighboring 
potential wells (see Fig. 1), with the states |0) and |1) 
of the qubit defined by the particle being in the left or 
right well (i.e. dual-rail encoding). A single quantum 
particle can occupy the two sites in a superposition, en¬ 
coding a qubit without the need for additional degrees 
of freedom. In this way, a system of n qubits can be 
realized in one dimension with n bosons and 2 n lattice 
sites, with one boson in the first two sites (representing 
the first qubit), one boson in the next two sites (repre¬ 
senting the second qubit), and so forth. Note that in this 
geometry, many physically permitted lattice states (e.g. 
those with more than one particle on the same site) are 
not members of the logical space (i.e. the multi-qubit 
tensor-product space). Nevertheless, we show that it is 
possible to engineer the lattice parameters such that, at 
time t = 1, U = e~ lH maps logical states only to other 
logical states with high fidelity, even though states out¬ 
side this subbasis are allowed at intermediate times. 


Implementing quantum gates — Having defined our 
qubits, we turn to the task of designing a universal set of 
quantum gates, i.e. finding lattice parameters that yield 
desired unitary transformations on the logical space. De¬ 
signing and building quantum logic gates remains one of 
the most difficult aspects of quantum computing, and our 
case is no exception. From the physical description of a 
given device — in our case, the lattice parameters — it is 
straightforward to write down the many-particle Hamil¬ 
tonian and to calculate the unitary evolution operator 



Figure 2: (color online). An implementation of the controlled- 
NOT (CNOT) gate according to the recipe in Eq. [5] (a) The 
real part and (b) the imaginary part of the two-particle uni¬ 
tary transform, U. The CNOT gate operation corresponds to 
the sub-matrix of the logic states, shown in solid-color bars 
and marked with red axis labels. Plots (c)-(f) show the posi¬ 
tion (in terms of the lattice sites, 1-4) of the particle density 
as a function of time, t, revealing the operation principle of 
the gate on each logical state (|00), 110), 101), and 111) respec¬ 
tively). One observes that the target qubit (in sites 3 & 4) 
performs Rabi-oscillations that are perturbed by the state of 
the control qubit (in sites 1 & 2) — the target qubit performs 
one fewer Rabi-flip if the control qubit is in the |1) state. 


U = e~ lH that fully describes the operation of the device. 
The inverse problem, however, is hard: given a desired 
unitary U, it is difficult to find a corresponding Hamilto¬ 
nian that meets the physical and geometrical constraints 
of the device, e.g. the one-dimensionality of the lattice. 
Furthermore, if the logical quantum states are only a 
subset of the full Hilbert space, then the quantum gate 
operation is only a sub-matrix of the overall evolution 
operator U. In this case, U is not even uniquely defined 
by the desired gate operation. As described below, we 
tackle these difficulties with a computational approach 
that finds appropriate lattice parameters to approximate 
a given gate operation with high fidelity. 

There are many options for the choice of a universal 
set of gates. One useful choice is the gate set of the 
controlled-NOT (CNOT) operation, along with either all 
single-qubit rotations (exactly universal) or the phase- 
shift gates and the Hadamard gate (approximately uni¬ 
versal) (21). In the following, we elaborate on the con¬ 
struction of the CNOT gate. Single-qubit gates involve 
only a single particle (i.e. no interaction terms) and are 
straightforward to calculate, as we detail in the Supple¬ 
mentary Information. 

To design the CNOT gate (a two-qubit gate) with the 
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dual-rail encoding, we consider a lattice with four sites 
and two bosons. This problem then is defined by eight 
lattice parameters: four on-site potential terms ( E m in 
Eq. 0, three tunneling terms and the interac¬ 

tion parameter (r). The complete two-body Hamilto¬ 
nian H is described by a 10 x 10 matrix (the size of the 
Hilbert space for two bosons in four modes). To perform 
the logical gate operation, the system is evolved accord¬ 
ing to U = e~ lH . The CNOT gate operation is then 
given by a 4 x 4 sub-matrix of U over the logical states 
11010) , 11001), |0110) , |0101) (presented here in the oc¬ 
cupation number basis); the other six basis states, while 
physically allowed, are not members of the logical basis. 
As explained above, finding the physical lattice parame¬ 
ters from the desired gate is a non-trivial inverse prob¬ 
lem. Using non-linear optimization techniques |25H28| 
(see Supplementary Information), we optimized the eight 
parameters of the system to maximize the fidelity of the 
gate when acting on the logical input states under the 
constraints that the parameters represent a physical one¬ 
dimensional lattice, i.e. that the on-site parameters are 
real, that the tunneling parameters are real and non¬ 
positive and connect only nearest-neighboring sites, and 
that the values of the on-site, tunneling, and interac¬ 
tion terms are within experimentally relevant bounds. 
Specifically, we demanded that 0 > Ji tTn > —J mQX , 
Jmax ^ ^ Jmaxi and T F max , where Jmax and 

r max are the largest allowed tunneling rate and interac¬ 
tion level in the optimization protocol. In our optimiza¬ 
tion we set J max = 4-7T, limiting the maximal number of 
tunneling events (or Rabi-oscillations) to 4. In practice, 
this experimental bound is dictated by the loss and de¬ 
coherence rate of the system, determining the maximal 
relevant propagation time. We also set T max = 10 J max . 
An example of a resulting lattice that yields the two-qubit 
CNOT gate is given (to two decimal places) by 


Gcnot = 7T 
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with interaction strength T = 21.687T. Here, the diagonal 
and off-diagonal entries of Gcnot represent the param¬ 
eters E m and J; jm , respectively, of the Hamiltonian H. 
Eq. [^represents a recipe for a four-site lattice that yields 
a CNOT gate with fidelity of 99.6%, whose operation is 
summarized in Fig. [2] 

If the bounds on the parameters are relaxed, the fidelity 
approaches even closer to unity. Fig. [3] summarizes the 
optimization results. Fig. |3](a) shows the convergence 
of independent runs with random starting points to the 
same final result. Fig. |3jb) presents the expected gate fi¬ 
delity vs. the maximally allowed values of the interaction 
r max . For a fixed maximal tunneling of Jmax = 47r, the 



Optimization step 


Figure 3: (color online). Optimization of the quantum-walk- 
based CNOT gate fidelity, (a) Convergence of different op¬ 
timization runs to the optimal gate fidelity. The solid black 
line and the shaded area represents the average and the stan¬ 
dard deviation values over 512 runs. Seven example runs are 
shown in the background (dotted lines), (b) Gate fidelity 
versus the maximum allowed interaction level F max , at a con¬ 
stant Jmax = 4-7T. (c) Gate fidelity for different maximal tun¬ 
neling rates Jmax at a constant maximal interaction level of 
Fmax — 207T. 


fidelity achieves a value close to 0.95 at T max /J max = 0.5 
and then slowly approaches unity as this value is further 
increased. In a system with a given r max , it is still possi¬ 
ble to improve the fidelity further by allowing more tun¬ 
neling events to take place, i.e. increasing J max ', see Fig 

HR 

Compiling a three-qubit primitive — Implementing a 
quantum algorithm using the scheme presented in this 
paper will involve several lattice configurations operating 
in sequence, as gates are sequentially applied in the algo¬ 
rithm. In principle, because the gate set presented in this 
work is universal, any multi-qubit operation can be bro¬ 
ken down into a sequence of single- and two-qubit gates, 
and thus implemented using the gates already presented. 
However, compiling common multi-step operations into 
a single primitive based on a single, time-independent 
Hamiltonian could reduce the possibility of errors aris¬ 
ing from dynamic changes to the lattice. As an example, 
we constructed a 3-qubit gate, shown in Fig. This 
gate is useful, for instance, in the 2-bit Deutsch-Jozsa 
algorithm |25], performing the oracle for the function 
f(x, y) = x ® y. (All other oracles for the 2-bit Deutsch- 
Jozsa algorithm are either a simple variation of this or¬ 
acle or require only single-qubit gates plus at most one 
CNOT gate.) Our computational approach allowed us to 
find a set of lattice parameters that realizes the complete 
three-qubit operation in a single gate. Fig. [4] presents 
an implementation of this three-qubit operation, at a fi- 
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(b) Imaginary (U) 

1 , 



frequency tu a ,m is detuned from the (now constant) cav¬ 
ity resonance w c by a frequency A m = oj a ,m — w c - The 
Hamiltonian is then well described by [40] 



(l,m) 



+ 9rn{al n a m + a+a m ), (5) 

m 


Figure 4: (color online). A 3-qubit operation of 2 CNOT 
gates (inset), compiled into a single gate U. (a) The real and 
(b) imaginary parts of our implementation of U. Only the 
logical basis states are shown. 

delity of 99.8%, using a single, one-dimensional six-site 
lattice: 
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( 3 ) 

with interaction strength T = 

108.24tt. 

In this case too, 


the fidelity could be improved by allowing larger tunnel¬ 
ing rates. 

Photonic systems — While we have focused our attention 
on cold atoms, it is possible to extend our scheme to 
nonlinear photonic systems, for example superconducting 
circuit QED systems m or nonlinear quantum-optical 
systems EM3- First, consider a system described by a 
similar Hamiltonian to that of Eq. 1: a set of coupled 
resonators exhibiting a single-photon Kerr nonlinearity 

m 

H = 'y ( cj m (o’h^am + — J 

m ^ ' 

+ 'y \ &m T 'y ~ n m (n m 1), (4) 

(l,m) rn 

where w m is the resonant frequency of the mth cavity 
and K the strength of the Kerr nonlinearity. Replacing 
Lo m with E m and K with T, we regain Eq. [I] exactly 
(neglecting the zero-point offset). Such systems have 
been experimentally demonstrated [ 35 ] an d site-by-site 
tunable coupled-cavity systems have been proposed [39] . 
However, our scheme requires a different E m (here, cav¬ 
ity frequency) at each lattice site, which is difficult to 
reconcile with a high-Q coupled resonant system. 

We can remedy this issue by instead coupling each of 
the cavities to a virtual two-level atom, whose transition 


where g m is the atom-photon coupling rate (i.e. vacuum 
Rabi frequency), and the a l m are the standard Pauli oper¬ 
ators. To second order in g m / A m , the state of the atom 
induces an a.c. Stark shift of the cavity by Xm = g„/A m 
m- Equivalently, a flip of the atom’s state shifts the res¬ 
onance of the cavity by 2\ m , meaning that any photons 
remaining in the cavity will have acquired an extra en¬ 
ergy of 2 Xm- Averaged over many Rabi-oscillations, this 
shifts the cavity photon frequency by Xm when two pho¬ 
tons are on the same site, causing a faster phase evolution 
rate than for a single photon. This is equivalent to the 
effect of an energy cost for two photons on the same site. 
Thus, by tuning g m and co m , we have enough degrees of 
freedom to obtain analogous dynamics to those of Eq. |T] 
having both distinct on-site energies and a tunable non¬ 
linear interaction. 

Of particular note, the inclusion of a two-level atom pro¬ 
vides a powerful advantage over both cold atom systems 
and the aforementioned nonlinear-Kerr cavities. As has 
been experimentally demonstrated m, due to the com¬ 
mutation of the a.c. Stark shift with the photonic and 
atomic states, it is possible to perform a quantum non¬ 
demolition measurement on the state of the atom (pho¬ 
tons) on a given site without affecting the photons (atom) 
on that site or any other states in the system. Potentially, 
this could be used to monitor or tune the operation of 
gates in situ or as part of an error correction scheme. 

Finally, we note that while we have focused on the ca¬ 
pabilities of superconducting systems in this section, de¬ 
vices based on quantum dots or graphene may offer sim¬ 
ilar opportunities at optical frequencies EIHST]. 

Conclusions — We have shown how quantum logic gates 
can be realized with high fidelity using the quantum walk 
of ultra-cold atoms on a one-dimensional lattice under 
experimentally achievable constraints. In particular, we 
gave a design for a high fidelity CNOT gate along with 
exact descriptions of single-qubit rotations, a computa¬ 
tionally complete set. Additionally, we demonstrated 
the compilation of a higher-order gate operation into a 
single operation. Our approach carries several impor¬ 
tant advantages over previous schemes. First, due to 
the dual-rail encoding we employ, the states of the sys¬ 
tem can be prepared and measured by simply placing 
and detecting single atoms at certain positions, both of 
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which are straightforward in present experimental sys¬ 
tems HU. Second, each quantum operation is carried 
out by a single, one-dimensional, time-invariant lattice 
potential. Third, the devices we propose are compact lat¬ 
tices of size 2 n (where n is the number of qubits) that can 
be realized on a line of potential wells with only nearest- 
neighbor hopping, in agreement with experimental capa¬ 
bilities. Our analysis shows that similar effects could be 
achieved in certain nonlinear quantum-optical systems. 
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Supplementary Information: Quantum Logic with Interacting Bosons in ID 


Single-qubit gates - To complement the presentation of a controlled-not (CNOT) gate in the main paper, we present 
the exact construction for a set of single-qubit gates that, together with the CNOT gate, make a universal quantum 
gate set. Since these are one-qubit operations, they are implemented using one particle in two lattice sites. As such, 
the interaction (T) is irrelevant and the matrix of lattice parameters 


( E x J 12 \ 

\ Jl2 E 2 ) 


(with J 12 < 0) can be directly interpreted as the Hamiltonian governing the single particle. The unitary gate, obtained 
by evolving with G for a time T, is then U = e~ lGT . 

One simple universal quantum gate set includes the Hadamard gate, the phase-shift gate, and the controlled-NOT 
(CNOT) gate (ST. The phase-shift gate is the simplest to implement. It is composed of two decoupled lattice sites in 
which the on-site energy between the sites is detuned. Specifically, to implement the single-qubit phase-shift operator 

Rg = f J j one may apply the single-qubit Hamiltonian 


Gr s = 




(SI) 


for time T = 1. (Obviously a different value for T may be chosen, so long as TGr 0 remains the same.) Next 
in complexity is the single-qubit Hadamard gate H = ^ j ^ 

generates the Hadamard transformation can be solved analytically, and is given by 


The Hamiltonian and propagation time that 


G h = 


V 2-1 -1 

-1 y/2 + 1 


T = 


2 V 2 ' 


(S2) 


Note that a simple tunnelling between two identical wells for half the tunnelling time (analogous to the beamsplitter 
operation in linear optics) will not reproduce the Hadamard gate in this case. Unlike the bulk quantum optics 
situation, in which the optical beamsplitter can be asymmetric in phase and thus can reproduce the Hadamard gate 
exactly, under our Hamiltonian dynamics the splitting is symmetric in phase and therefore modified tunnelling rates 
and additional diagonal terms are required in order to correct the output phases. We note this holds true also for 
integrated quantum photonics gates jS2]. A summary of the physical construction for the single-qubit operations is 
given in Fig. [Si] 

Alternatively, one can use the (exactly) universal gate set of CNOT together with all single-qubit unitaries. Any 
single-qubit unitary, U, can be implemented by first decomposing it in the form (Slj 


where the 2 -rotation R z (9 ) 


U = e ia R z (p)R x (j)R z (6) = e ia R z (fi)HR z (ry)HR z (5) 
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and the a;-rotation R x {9) = exp 
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can be implemented either with the Hamiltonian 
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(S3) 


or by conjugating R z (6 ) by the Hadamard operation H described earlier. The phase e*“ can be implemented with 

G a = -a(jj),T=l. 

Quantum process tomography for CNOT - The basic idea of quantum process tomography jS3. is to completely 
determine the mathematical structure of a given operation (in our case, the CNOT on two qubits), by applying the 
operation to a variety of different input states and measuring each of the resulting states in a variety of bases. 
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Figure SI: (color online). A double potential well containing a single-boson configuration representing a single qubit in the (a) 
|0) state and the (b) |1) state, (c) An implementation of the single-qubit phase gate, according to the recipe in Eq. SI (d) An 
implementation of the Hadamard gate, according to the recipe in Eq. |S2| 


Tomography procedures allow much freedom is choosing input states and measurement bases. We want to make such 
choices wisely to facilitate performing the procedure on our particular experimental apparatus. Two main points guide 
these choices. One is that the only two-qubit gate be the purported CNOT gate itself — all other gates required for 
the tomography should be relatively simple so that we adequately assess the CNOT as the major source of errors. 
The second is that the other gates be easy for the Bose-Hubbard model to implement, preferably with as few gates 
as possible (since limitations in coherency times may limit the number of gates reasonably performed using current 
technology). Here we suggest input states and measurement bases to accomplish these goals, requiring only two single¬ 
qubit Bose-Hubbard Hamiltonians (one prior to the CNOT, and one following it) per (computational) measurement. 
Specifically, each measurement in the tomography involves applying an operation of the form (V 3 OV 4 ) CNOT (ViOV^) 
to a computational basis state, followed by a measurement in the computational basis, where each V) can be performed 
with a single one-qubit Bose-Hubbard Hamiltonian; we therefore need only keep the state coherent for three consecutive 
operations. 

We assume that the reader is following the quantum process tomography procedure outlined in Ref. [Si] , which requires 
measuring Pauli observables. To measure a Pauli matrix a € {/, X, Y, Z}, one must perform a basis-transformation 
gate to their output state so as to measure in the eigenbasis of a. For a = I and Z, this is trivial, as the computational 
basis suffices; however, for a = X and Y, we recommend measuring in the following eigenbases, which can be performed 
by first applying the following basis-transformation matrices, 


and 
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(S4) 


(S5) 


Observe that Ux is the Hadamard matrix, the Hamiltonian for which was given in Eq. S2 and that Uy = R x (ir/2), 
which is generated by the Hamiltonian in Eq. S3 with 9 = 7 r/ 2 . Note that for Y we chose a non-standard orthonormal 
eigenbasis so that the basis-transformation can be accomplished using a single Bose-Hubbard gate . 


The input states required to perform the quantum process tomography procedure of Ref. |S1| are straightforward to 
produce. For example, one may use tensor products of the following single-qubit input states: the computational basis 

states | 0 ) and 1 1 ), the +1 X-eigenvector |+) = Ux | 0 ) = -j= f j j, and the +1 T-eigenvector |—) = Uy |0) = 
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where Ux and Uy are given in Eqs. (S4| and (S5). Note that Uy 
in Eq. S3 with 8 = 7n/2. 


Rx{7Tt/2) can be generated by the Hamiltonian 


Computational Methods - In this section, we detail the numerical methods used to find the gates presented in the main 
paper. We made extensive use of a free-software implementation of a variety of numerical optimization algorithms fS4| . 
This enabled us to, with a single specification of cost function and constraints, compare the success and computational 
cost of a number of different optimization approaches. We discovered that a randomly-seeded global optimization 
algorithm (S5l, IS6| combined with a gradient-free local algorithm |S7] gave the best performance, both in terms of 
number of iterations and computational run-time. 


Careful selection of the cost function was crucial to the success of this work, and interacted with the choice of the 
aforementioned algorithms, particularly the local optimizer. Throughout the paper, we define the fidelity of the gate 
in terms of the Hilbert-Schmidt inner-product between the target unitary gate operation Uq and the unitary operation 
U generated by the Hamiltonian at a given step of the optimization (restricted to the logical subspace). Specifically, 
the fidelity is defined to be 


F(Uo,U) = \{Uo,U) c \ 


with 


(Uo,U) c = 


Tr(U\U) 

N 


where N is the dimension of the logical space (4 for two-qubit gates). This fidelity can be interpreted as a sort of 
lower-bound average fidelity of the gate. 

To be specific about the iterative numerical process, at each step we generated U from a vector corresponding to 
lattice parameters and calculated F(Uq, U ). Numerically, we found that minimizing the function 1 — F 2 , rather than 
1 — F, gave superior performance. In the case of the algorithm given in Ref. |S7j . the reason for this is clear: the 
algorithm assumes a quadratic cost function. However, we found that even with algorithms designed for linear cost 
functions (e.g. [[S8j). convergence was much slower than for the quadratic cost function. 

Finally, in order to ensure that U has the same global phase as Uq (this is for aesthetic purposes, as F(Uq, U) is 
invariant under multiplication by a global phase), we placed a cost on the phase of the matrix element u\p. We 
found this to be most efficiently implemented by adding the term sin(arg(rti : i)) 2 to the cost function. This function 
is quadratic when perturbed about zero, is non-negative, and is symmetric about nn for all n £ Z, making it an 
ideal candidate function. We verified that the introduction of this additional cost both yielded a U with appropriate 
phase (see Fig. 2 in the main text) and did not result in a decreased fidelity compared to optimization without this 
constraint. 
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